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APPLICATION  OF  A  FEATURE  SELECTION  TECHNIQUE 
TO  SAMPLES  OF  HIGH  RESOLUTION  SYNTHETIC 
APERTURE  RADAR  IMAGERY 


The  purpose  of  this  research  note  is  to  show  the  appli- 
INTRODUCTION  cation  of  a  feature  selection  technique  to  samples  of 
synthetic  aperture  radar  imagery  and  to  present  some 
preliminary  results.  In  the  past,  feature  selection  techniques  have  been  applied  to 
data  computed  from  digitized  aerial  photography.  However,  it  appears  that  no  one 
has  as  yet  applied  feature  selection  techniques  to  high  resolution  radar  imagery. 
In  order  to  perform  classification  of  terrain  features  using  radar  imagery,  feature 
selection  is  an  important  initial  step.  Feature  selection  consists  of  choosing  those 
features  that  are  most  effective  for  showing  class  separability  and  for  performing 
a  reduction  in  the  dimensionality  of  the  feature  vector.  The  following  sections  will 
present  a  discussion  of  the  feature  selection  technique,  along  with  its  application 
to  selected  samples  of  radar  imagery. 


The  application  of  pattern  recognition  techniques  is 
METHODOLOGY  accomplished  usually  in  two  steps,  namely,  feature  selec¬ 
tion  and  classifier  design.  The  feature  selection  process 
that  precedes  the  classification  process  consists  of  techniques  applicable  to  one 
class  or  to  multiple  classes.  The  feature  selection  technique  to  be  discussed  here 
is  applicable  to  multiple  classes.  It  provides  the  capability  of  reducing  the  number 
of  components  of  the  original  feature  vector  in  such  a  way  that  the  resulting  com¬ 
ponents  are  optimized  to  show  class  separability.  The  feature  selection  technique 
comes  from  the  field  of  discriminant  analysis  of  statistics  and  is  independent  of  the 
probability  density  functions  of  the  feature  vector  data.  The  feature  selection 
operation  can  be  expressed  as  a  linear  transformation  of  the  following  form. 

Y  =  AX  (1) 
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where  X  is  the  original  feature  vector  with  dimensionality  n  x  1 ;  A  is  the  trans¬ 
formation  matrix  of  dimensionality  m  x  n,  where  m  is  less  than  n;  and  Y  is  the 
transformed  feature  vector  with  dimensionality  m  x  1.  The  feature  selection 
problem  is  now  reduced  to  determining  the  matrix  A.  In  order  to  calculate  A, 
use  is  made  of  the  within-class  and  between-class  scatter  matrices.  A  within-class 
scatter  matrix  shows  the  scatter  of  samples  around  their  class  expected  vector  and 
can  be  expressed  by 


N 

sw  =  £  PfcOjiCj  (2) 

i=  l 

where  Sw  is  the  within-class  scatter  matrix.  P(cOj)  is  the  a  priori  probability  of 
the  ilh  class,  Ct  is  the  covariance  matrix  of  the  ith  class,  and  N  is  the  total 
number  of  classes.  A  between-class  scatter  matrix  can  be  defined  in  many  ways; 
however,  the  following  definition  was  the  one  utilized  here: 

Sb  =  C,  +C2 +(M, -M2)  (Mj -M2)t  (3) 

where  Sb  is  the  between-class  scatter  matrix,  C,  is  the  covariance  matrix  for 
class  1,  C2  is  the  covariance  matrix  for  class  2,  M,  is  the  mean  vector  for  class  1, 
M2  is  the  mean  vector  for  class  2,  and  T  means  transpose.  This  definition  of  the 
between-class  scatter  matrix  is  valid  only  for  the  case  where  N  (the  number  of 
classes)  is  equal  to  two.  In  order  to  have  criteria  for  class  separability,  a  number 
must  be  derived  from  the  scatter  matrices.  This  number  should  increase  when  the 
distances  between  points  belonging  to  different  classes  are  increasing  or  when  the 
distances  between  points  belonging  to  the  same  class  are  decreasing.  One  criterion 
is  the  use  of  J,,  which  can  be  defined  as  follows. 

J.  =  trace  (S  1  S  )  (4) 

*  2  m  1  m 

where 


S,  =  AS  A1  and  S.  =  AShAT 

2  m  w  i  m  u 
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The  feature  selection  problem  now  requires  that  we  select  the  particular  trans¬ 
formation  matrix  A,  which  maximizes  the  value  of  Jj.  Fukunaga1  shows  that 
A  is  made  up  of  the  normalized  eigenvectors  of  the  matrix  Sw*  Sb. 


At  =  [0,  <t>2 . <t>m  1  (5) 

where  4>x  is  the  eigenvector  associated  with  the  largest  eigenvalue,  <t>2  is  the 
eigenvector  associated  with  the  second  largest  eigenvalue,  etc.  Once  the  matrix  A 
is  computed  from  (5),  the  new  feature  vector  Y  can  be  computed  for  each  point 
in  each  class. 


The  feature  selection  technique  was  applied  to  samples 
INVESTIGATION  of  high  resolution  synthetic  aperture  radar  imagery  taken 
over  the  Huntsville,  Alabama,  area  with  the  APD-10 
radar  system.  Sections  of  the  radar  imagery  were  digitized  and  stored  on  a  digital 
disk  unit.  A  Lexidata  system  3400  display  processor  was  used  to  display  the  images 
on  a  cathode  ray  tube  and  to  take  100  samples  for  each  of  four  terrain  classes 
from  the  imagery.  Each  sample  consisted  of  a  32  by  32  pixel  element  window 
located  within  a  section  of  one  particular  terrain  class.  The  four  classes  considered 
were  cities  (combination  of  commercial  and  residential  structures,  DLMS  cate¬ 
gory  $504  FIC  301  and  $505  FIC401),  fields  (agriculture  used  primarily  for 
crop  and  pasture  land,  DLMS  category  $5 1 0  FIC  950),  water  (rivers  with  smooth 
fresh  water,  DLMS  category  $510  FIC  940)  and  fresh  water  subject  to  ice  (lakes 
and  reservoirs,  DLMS  category  $510  FIC  943),  and  forests  (mixed  trees,  deci¬ 
duous  and  evergreens,  DLMS  category  $510  FIC  954).  A  feature  vector  con¬ 
sisting  of  13  components  was  computed  for  each  sample.  These  13  components 
were  made  up  of  the  first-  and  second-order  gray  level  histogram  statistics  computed 
from  each  sample  window.  The  explicit  equations  used  for  the  13  components  of 
the  original  feature  vector  are  shown  in  appendix  A.  A  computer  program  was 
written  for  the  Hewlett-Packard  1000  computer  to  calculate  J,  as  a  function  of 
displacement  in  x  and  y.  This  computer  program  was  also  used  to  calculate  the 
transformation  matrix  A.  A  listing  of  this  program  is  provided  in  appendix  B. 
A  second  computer  program  was  written  to  calculate  the  new  feature  vector  Y  for 
the  four  hundred  samples  taken  from  the  radar  imagery.  A  listing  for  this  second 
computer  program  is  given  in  appendix  C. 


*Keinosukc  I  ukunaga.  Introduction  to  Statistical  Pattern  Recognition.  Academic  Press,  1972. 


In  this  section  some  results  of  numerical  calculations  are  pre- 
RESULTS  sented.  Table  1  shows  the  results  of  calculating  the  value  of 
J  j  for  various  values  of  displacement  DX  and  DY  for  the  two 
classes  of  forests  and  cities.  The  a  priori  probabilities  for  the  two  classes  werfe 
assumed  equal  to  O.5.2  In  table  1  the  largest  value  of  Jj  is  29.8525,  which  is 
associated  with  a  DX  of  -3  and  a  DY  of  4.  The  significance  of  a  maximum  value 
of  lx  occurring  at  these  particular  values  of  displacement  is  not  understood  at  this 
time.  The  values  of  the  displacement  DX  and  DY  associated  with  the  largest  value 
of  Jj  were  utilized  to  compute  a  new  feature  vector  Y  with  two  dimensions. 

Figure  1  shows  a  plot  of  two  second-order  statistical  components  for  the  forest  and 
city  samples.  The  triangles  are  the  results  from  city  calculations;  and  t'  x’s  are  the 
results  from  forest  calculations.  This  figure  clearly  shows  the  neces*  tor  feature 
selection  because  the  original  data  for  forests  and  cities  is  not  sepai  \  Figure  2 
presents  the  results  after  the  transformation  Y  =  AX  is  applied  vie  original 
feature  vector  X.  In  this  figure  the  x’s  again  represent  calculatio’  >m  forests. 
We  now  see  that  the  data  from  forests  is  clustered  and  almost  totally  ed  from 

the  data  for  cities,  represented  by  the  triangles.  Even  though  a  jints  still 
overlap,  the  improvement  is  very  dramatic.  The  transformation  has  succeeded  in 
reducing  the  dimensionality  of  the  original  feature  vector  from  1 3  to  2  and  also  in 
separating  the  clusters  of  data  for  the  two  classes. 

Figure  3  shows  the  results  of  the  transformation  for  the  two  classes  of  cities 
and  fields.  The  x’s  represent  the  data  from  cities.  Good  separation  is  obtained  for 
the  two  classes  as  only  a  few  points  are  overlapping.  Figure  4  shows  the  results  of 
calculations  for  the  two  classes  of  cities  and  water.  The  x’s  represent  the  data  from 
cities.  In  this  particular  case  the  two  classes  are  totally  separated  as  well  as  being 
clustered  fairly  well.  Figure  5  presents  the  results  for  forests  and  water.  The  x’s 
represent  the  data  for  forests.  Again  the  two  classes  have  clearly  been  separated  by 
the  transformation.  Figure  6  shows  the  results  for  fields  and  water,  with  the  data 
for  fields  represented  by  the  x’s.  These  two  classes  have  also  been  totally  separated. 
Figure  7  shows  the  results  for  forests  and  fields  with  the  data  for  forests  repre¬ 
sented  by  the  x’s.  The  two  classes  are  not  well  separated  and  other  features  may 
have  to  be  investigated  to  obtain  better  separation.  Another  possibility  would  be  to 
try  using  nonlinear  feature  selection  techniques  for  these  two  classes. 


^Thit  assumption  of  equal  probability  for  cities  and  forests  appeared  to  be  appropriate  for  the  area  where  the  samples 
were  taken. 
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TABLE  1.  Values  of  J,  for  first  and  second  order  histogram  statistics 


DX 

DY 

J, 

1 

0 

22.1485 

2 

0 

19.8078 

3 

0 

19.4333 

4 

0 

19.3099 

5 

0 

19.7508 

6 

0 

20.3677 

7 

0 

19.8529 

0 

1 

21.0314 

0 

2 

21.76-3 

0 

3 

21.3524 

0 

4 

21.7343 

0 

5 

22.3144 

0 

6 

23.1462 

0 

7 

23.5903 

1 

1 

20.6771 

1 

20.4869 

2 

24.3979 

3 

2 

23.9334 

4 

2 

21.5577 

3 

3 

25.7921 

-  1 

1 

20.2745 

2 

24.4373 

-  3 

25.0272 

-  3 

4 

29.8525 

4 

4 

22.7131 

1 .  The  feature  selection  technique  discussed  in  this  report 
CONCLUSIONS  appears  to  be  a  powerful  tool  for  the  application  of 
pattern  recognition  to  high  resolution  synthetic  aperture 
radar  imagery. 

2.  In  order  to  separate  forests  and  agricultural  fields  it  appears  that  the  first  and 
second  order  histogram  statistics  combined  with  a  linear  feature  selection  technique 
may  not  be  sufficient  as  a  feature  vector. 


FIGURE  7.  Forests  and  fields  DX 
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APPENDIX  A 

FEATURE  VECTOR  COMPONENTS 

The  following  first  and  second  order  histogram  measures  were  used  to  construct 
a  thirteen  dimensional  feature  vector: 


L— 1 

Mean 

b 

-  E 

bP(b)  =  x  j 

b  =  o 

Variance 

"2b 

L  —1 

-  E 

(b — b)2  P(b)  =  x2 

b  =  o 

i 

L—  1 

Skewness 

bs 

o 

II 

£  (h-h)3  p<b>  = 

b  =  o 

Kurtosis 

bK 

1 

L-l  _ 

v  a>-b)4  P(b) 

b  =  o 

Energy 

bN 

L  —1 

=  E 

b=  o 

(P(b)| 2  =  xs 

Entropy 

b. 

II 

P(b)log2  [P(b)l  =  x6 

b  =  o 

Autocorrelation 

Ba 

L-l 

=  E 

L-l 

y-'  abP(a.b)  =  x7 

a  -  o 

b  =  o 
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L-l 

L-l 

Covariance 

Bc 

=  E 

E 

(a-a)  (b-b)  P  (a,b)  =  xg 

a-o 

b  =  o 

L-l 

L-l 

Inertia 

B, 

=  E 

E 

(a-b)2  P  (a,b)  =  x9 

a  =  o 

b  =  o 

L-l 

L— 1 

Absolute  Value 

By 

=  E 

E 

|  a-b  |  P  (a,b)  =  x10 

a=  o 

a  =  o 

L-l 

L— 1 

Inverse  Difference 

B» 

=  y 

y 

P(a.b) 

Ld 

a  =  o 

Ld 

b  =  o 

1  + (a-b)2  Xl1 

L-l 

L-l 

[P(a.b)l 2  =  x j  2 

Knergy 

Bn 

=  E 

E 

a-o 

b  =  o 

L-l 

L-l 

Lntropy 

B„ 

=  -E 

E 

P(a.b)log2  |P(a.b)|  =  x,3 

a  '  o 

b  =  o 

whore  L  is  the  number  of  grey  levels  and  P(b)  and  P(a,b)  are  given  below 


P(b) 


Q(b) 

M 


M  is  the  total  number  of  pixels  in  the  sample  window.  In  this  ease  M  was  equal  to 
1024.  Q(b)  is  the  number  of  pixels  of  grey  tone  b  which  occur  in  the  sample  window. 

O(a.b) 

Pla.b)  =  ___ 

Q(a.b)  is  the  number  of  times  greytone  a  is  located  next  to  greytone  b  by  the  dis¬ 
placement  Ax  and  Ay. 
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APPENDIX  B. 

Computer  Program  for  Calculating  J,  and  the  Transformation  Matrix  A 
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1)  3  1  S 

h  J  thM 

03U* 

00  <!0  K !  NOWAY 

»  3  l  v 

;jo 

r:tiM  t  <k  >  -  o .  o 

0  3  .  M 

00  30  1  -- 1  3t.‘ 

II 3 »  v 

DO  30  J • 1 . 3? 

0  3,  mi 

I  K  -  f  AW  <  J  J  ) 

O.V  i 

fJON  M  1 W  +  l  i  TON  1 

0  5,  V 

30 

UiN  I  j  NOl 

0  3/ » 3 

Xh  -  M 

.i  V* 

oo  ys  k  -  »  .  noway 

n  vji. 

no  m  i  noway 

0  V.’t, 

* : (M.I  l  <  K  .  M  >  0  .  u 

0  S3  • 

MX  .3^3  1  OX 

l»3/H 

1  F  i  3.? -NX  )  !  r  .14 

0  3  ,J  V 

1  /? 

J.  NOW  -MX-.V 

03  ?n 

NX  0-1 

0  5  5  i 

NXF  -3,-f  -  )  N(  lv 

0  3  a, • 

«  0  10  1 v 

11333 

1  4 

NX  H  -  1 

03  4  *4 

NX1_-  V 

1)3  5^ 

i;o  hi 

0  *,  *)•% 

1  /.» 

i  Ni  W  5 /-'-Pi x 

0  33/ 

N  X  H  1  1  +  1  N/  W 

0  3  %D 

iUl  -  „s3  NXO 

033v 

IV 

M  r  3iJ  ♦  1  />  * 

ii  vHI 

ll  i  JtJ  Ml  i.'l  «  3  . 

034  i 

1 

i  i'll  W  -  MY  3.- 

o  5  * 

i'1 1  |i-  i 

II  54  s 

Nil'  3  3  1  Of  1. 

ii  *,  4 «) 

in  ni  . '  - 

II  5  4'-. 

J3 

N  v'F<  -  1 

ii  i 0 1 . 

Nil  -  3t‘ 

M34  *' 

Mi  Hi  • 

ii  *,  1 1 . 

.  '» 

i  Ni  R  n  > 

ll  A  <V» 

Mi  H  \  *■  1  N  f .  W 

II  .  U 

■ i  1  -  si  •  N »  | ; 

H  5  3  i 

I'lM  1  1  NUr 

•1  5  ■ 

Hll  \l,  |  |  IV  1  i 

Ii  *i  i 

Oil  \  >  1  i  .NO 

•  1  >l ,  •* 

1  ,  .  1  mK  .  It  ,  li  f  1  ) 

ll.i  . 

l  .  | AW ( 0 i H  )  I  1  ♦ 

M 

. .  t  1  i  «  .  1  .  ♦  i  ' 

II  *  '  i  . 

A  . 

i  nr'  i  i  M  h 

1  6 


N»)(*  .1  )  i 


-  •  I  Mi  I  l  I  1  :  ) 


i  ♦  *  :  ■  •  - 
i  . 


mii  . J  1 1  U  i  .VI.I.  HI 
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APPENDIX  B 
(Continued) 


li  3S9 

40 

P(K  >--=CONT<K  >/XM 

0  360 

0(1  4S  K=1  NGRAY 

0  36  l 

4S 

AVG-AVG+<K-t)*P<K> 

0  36,! 

1>U  SO  K-l .NGRAY 

0  363 

SO 

9AR-9AR  +  <  <K-1-AUG)#*?)*P<K> 

0364 

SI  D-sWNl  <UAR> 

0  J6S 

DO  SS  K-l • NGRAY 

0366 

ss 

SKU-SKW-MU  1-AUU1**3)*P(K> 

0  36  V 

8KW-'iKU7ST0**J 

0  366 

OU  60  K-l . NGRAY 

11369 

60 

XK  1  -XK  I ♦  (  <K-1-AUG)**4>*P(K) 

0  170 

XK  1  =■■<  XK  1  /  <s;i)**4)  )  -3  . 

03  v  i 

DO  6S  K-l  .NGRAY 

0  37.* 

6S 

1  Nl'.-l  NG‘P  <  K  I  *P  <K  ) 

0  37  3 

Dll  70  K  - 1  .NGRAY 

0  374 

1  p  ( P  (  K  >  1  (JO.  0  )  GO  TO  70 

037S 

KM  1  -EN  f  +  3. 331929*8  IK  >*AI  OCT  (  P  (  K  1  1 

0376 

70 

I.UNI  JNUE 

11377 

ENI --FN1 

0  37 ti 

|)(l  7S  K  “  1  .  NGRAY 

0  379 

1)1)  7S  8*1  NGRAY 

0  3HU 

7S 

PPIK  .  8  1  -1.011 1  IK  .  8  )  /XM 

0  381 

Dll  80  K-l  NGRAY 

0  382 

1HI  80  H-)  NGRAY 

0  383 

All  I  -AUIMK-1  >*<8-l  ItPPIK  .8) 

0  384 

80 

(UNI  1NII1 

0  18S 

AOr.K  -0,0 

0  386 

AVGN-0 . 0 

0387 

DO  HS  K  l . NGR A) 

0  388 

Dll  8S  M- 1  .NGRAY 

0  389 

AOCK - AOGK ♦  I  K  -  1 >*PPIK  .H) 

039  0 

8S 

AOGN- AVGM+  <8-1  I *PP  <  K . 8  I 

0  391 

t)(J  90  K-l.  NGR  AY 

03v;? 

DO  90  8-1. NGRAY 

0.393 

90 

X1F.  =  XU  +  <K  8>*<K-8)*PP<K -M> 

0  3V4 

DO  9S  K-l .NGRAY 

0  3  VS 

DO  9S  8-1  NGRAY 

0396 

APSS=A6SS*< jAPSIK-8)  >*PPIK  .8) 

0397 

9S 

CON  1 1  NOE 

0398 

DO  100  K - 1 .NGRAY 

0399 

DO  100  8* 1 .NGRAY 

0400 

100 

X  J  ll-Xi  D+PP  <K  .81/(1  .  O+IK-M )«*<!) 

040  t 

00  10S  K-l. NGRAY 

0  402 

DO  10S  8- 1. NGR AY 

0403 

10S 

E  NY  -ENY+PP  I  K  . 8I4PPIK  .8' 

0404 

DO  1 10  K-l .NGRAY 

040S 

DO  110  8=1 , NGRAY 

0406 

IEIPPIK  .8)  ,E0.  0  .  OIGO  (0  110 

0407 

FNR=KNR+3. 3219294PP <K . M 14ALOG I (PP ( K .8) ) 

0408 

110 

CONI J NOE 

0409 

ENR=-FNR 

0410 

DO  1  IS  K  =  1  .NGRAY 

0411 

DO  US  8=  1  .NGRAY 

0412 

i:i>V=l.(JV+  <  K  - 1  -  AVGK  )  *  <  8-  1  -  AVGH  >  *PP I K . 8  > 

0413 

US 

CONTINUE 

0414 

RETURN 

041 S 

FND 

041  6 
0417 

END* 

I 
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■WAX  I  R 

I)  0  0  4 
o  (i  n  2 
0003 
0004 
000b 
0006 
0  0  07 
0  0  08 
000V 
00  to 
001 1 
0  0  1  2 
0  0  V  3 
0014 
001b 
0016 
0017 
0  018 
0  0 IV 
0  020 
0  021 
0021? 

0  023 
0  0  24 
0  02b 
0  036 
0  027 
0028 
0  U2V 
0  0.10 
0  031 
0  032 
0  033 
0  0  34 
0  03b 
0  036 
0  037 
0  038 
0  037 
0040 
0  0  41 
0  0  42 
0043 
0044 
0  0  4b 
0046 
0  0  47 
0  048 
0  04V 
OObO 
OflSl 
0  0S2 
0  Ob3 
0WS4 
OObb 
OOS6 

0  osv 

0  0S8 


APPENDIX  C. 


Computer  Program  for  Calculating  the  Transformed  Feature  Vector  Y  =  AX 


T  “  00  0  0  4  IS  (JN  CRUOOli  US  I  NO  00013  HI  K  !■  R-0000 


F  I  N4  ,  L 


ixccccuxccixeccccccccccccccccccixcccccccccccccccccccucccxcccccccccccc 

C  (HIS  PROGRAM  COMPUTES  THE  TRANSFORMA I i ON 


0  Y-AX  WHERE  A  18  COMPUIEP  F  ROM  ANO I  HER 
C  PROGRAM  AND  IS  INPUT  HERE 

(Xcmxixcccixci.cccxixcxixxccccxixccxecccrcccccixccccccixcxixixLxcixrx 


PROGRAM  YAX IR (3. tOOO ) 

DIMENSION  (.001  ( 5 1  .  I  MAGE  <10241  . X < I  00  ,  13)  .  Y(  1  00 .2)  .  A<2 , 1 3) 
DIMENSION  (008(144) . IE  UE  (3  >  .  AM <  2 ) . OAR  < 3 ) 

DIMENSION  AV<  6.6)  , D  <  3 . 2  > 
tool VALENCE  < IMAGE . Y > 

CALL  RMPAR  <  LUO  I ) 


I  II  -‘LUO  I  <  1  ) 


CALL  (RLUdll) 

N02 

4  K1RMA T<  1 1  ) 

WRITE <LO, 6) 

6  FORMAT  <  “tNIt.R  I  HE  VALUE  FOR  M  I  HE  NUMBER  OF  ROWS  (N  I  HE  IRANSFOR 
1MAIION  MATRIX  A") 

RFAD<1  0.41M 
WR l IF  <LU ,8) 

8  FORMAT ( "ENTER  THE  TRANSFORMATION  MATRIX  A") 

DO  10  l  »=  1  .  M 
DU  10  J-i , 1 3 


WRITE  (HI.  12)1  ..I 

12  FORMAT < “1“" . 12.2X. "J-“ . J 2 , *A< I . J >«* ) 

READ<  LIJ  ■  1  4  )  A  ( 1  .  J  ) 

14  FORMA) <E V . 3) 

10  CONTINUE 
WRITE (6.7V) 

VV  T  URMAX  1  X  .  "  CHE  TRANSF  ORMATION  MATRIX  A“) 

DO  78  I - 1  .  M 

WR1  Tt.<6 .77)  <A< )  ,J)  ,J  =  1  ,13) 

77  FORMAT) IX . 1 3  <  E V . 3 . IX) ) 

78  CUNITNUE 

DO  80  NJ--1  .  NC 
WR  J  I  E  <  1  U  .  1  6  ) 

16  FORMAT  <  "ENTER  THE  NUMBER  OF  IMAGE  SAMPl  ES  TO  BE  ANAl  Y/ED  . IE. 100") 
Rt  AD  <  LU .18) NDATA 
18  FORMAT (13) 

UR  I TFU.U .20 ) 

20  FORMAT < "ENTER  THE  FILE  NAME  FOR  I  HE  DA  I A  SET") 

READ  (LU  .  2?  )  IE  II  F 
22  F  ORMA I <  3A2 ) 

WRITE (l.U.  IS) 

lb  FORMA  I  ("DISK  t.O  NUMBER-)") 

RE  AD  <  LU  .2101)  1  DLl) 

WRITE  (LIJ.  21  00) 

2100  FORMAT < "ENTER  A  VAl Ut  FOR  1DX”) 

READ  <  I..U  .210  V  )1DX 
3101  FORMAT (12) 

WR  T  IF  <  L.U  ,  230  0  ) 

220  0  FORMAT  <  "ENTER  A  VAt  UE  FOR  (DY  "  1 
READIED.  21  01  )IDY 
WR  I.  T  E  (  6 . 1  400  )  I  DX  ,  1  DY 
1400  FORMAT (2X . "IPX*" .12. SX, "IDY-" .121 
WRTTE:<6.200)  (KILE 
20  0  E  OR MAT (IX. 3A2 ) 
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APPENDIX  C. 
(Continued) 


0  0  5V 
00*0 
00M 
0  0*? 
0  06.4 
00*4 
0  0*5 
00*6 
0  0*7 
0  0*R 
0  0*9 
0070 
0071 
0  0  7  7 
0  07  A 
0  0  74 
0  075 
0  07* 
0  077 
0  0  78 
0  0  79 
0  080 
0081 
0  08? 
0  08-4 
0  084 
0  085 
0  08* 
0  087 
0  088 
0089 

o  o  v  n 
0  091 
0  0  9? 
0  096 
0094 
0  0  95 
0  09* 
0  097 
0098 
0  099 
0  1  0  II 
0  1  0  1 
(110? 
0  I  0  .4 

I)  )  II  4 
OKI  5 
I)  I  0  * 
0  l  0  7 
0  l  08 
II  1  0  9 
0110 
0  1  !  I 

II  1  I  7 

II I  i  1 
0  114 
HI  i'i 
0)1* 
0  I  i  7 
II  I  I  8 


CALL  OPF.Nl  1  DC  8 . 16KR  .IHU.0,11  ,-llH.U) 

IF  (  l'fc'RR  .LI.  OIGIJ  TIJ  7000 
GO  10  74 

7000  UR  I  TE  ( Lll .  201 1)  >  IFRR 

2010  F ORMAT <  "OPt  N  FJl.t  1  RRUR  "  .15) 

GO  TO  999 
24  )Cl)NT  =  l 

CAU  LA8TN( I  DC 8 . * ) 

1  .1  ,1-1 

DO  19  1:1.8 

CALL  RF ADI  ( 1  DC 8 . 1LRR .  I MAGL  <  J  > ) 

19  J-J+128 

1 1  <  1  F.RR  .11.  0  ICO  III  .1000 
Cl)  TO  26 

.100  0  I  CON  I -NDA I  A 

UR  ITt  <  Lll  -  2070  )  IFRR 
7020  KJRMA  I  (  "RFAD  F  ILF  IRROR'MS) 

CO  TO  999 
7*  NCR AY  =  16 

CAl  1.  lOCAl  <  1MACL  .  IMAGE  .  1  024 .0.15) 

CAL  1.  FFOFCI  I  MAC!  .  AVC  ,  OAR  .  8KU  .  XK  I  .1  NG  ,ENI  .  AO  I  . CUV  .  X  1C.  ,  AKSR  ,  X 1  D  . 
IF.NY.6NR  .  NCRAY  .  1  DX  .  1I)Y  ) 

XU  CONI'  ll-AOG 
X( I CON  1  . ? ) -VAR 
XI  (CONI  ,  .5  )  -  8KU 
X(  ICON  I  .4) -XK I 
x (  ICON  I  .  5  >  F  NC 
XI  I  CON  I  .*)  INI 
X< ICON  I  -Vl-rtOT 
X<  1  CON  I  .  8) --CIJO 
X( ICON  I  .9)  -  Xll 
X(  ICON  I  .  1  OK-AHSS 
x<  I  CON  1.11  )  *■■  X 1 1  > 

X< ICON  I  . t?>u NY 
X  C  ICON  I  1  .4  )  NR 
IF  (  ICON  I  -NDA  I  A  178  .  40  .  4  0 
70  (CONI  -  ICON l M 
1.0  10  IT 

-.1)  UR  J  ltd  0  .  -4 . '  I  N  ,1  ,  IF  II  l 
.47  FORMA  I  <  "N.l:-"  ,  11  ,.<X  .  ,<A7  > 

Dll  -44  j  -  i  NDA  I  A 
1)0  ,44  RUN 
.44  t  <  i  .  K  )  -  0  .  0 

1)0  .4*  r  I  .  NDA  I  A 
Dll  ,4*  K  !  M 
Dll  .4*  MK  -  I  .  t  ,4 

.4*  ill  K  )  ■  1  (  I  .  K  )  i  A(  K  ,MK  >*X<  1  ,  MK  ) 

CAl. I  CPI  III  (  Y(  l  .  1  )  .  Y<  1 .7)  .  1  IHi  III. Nil 
I'll  Ml  K-l  M 
AM  I  K  > "-  U  .  Ii 

nK-ndai a 

DO  45  .1  I  .  NDA  I  A 
4  5  AM  <  K > -AM  I K  ) ( 1  (j  ,K 1 
511  AMI  K  )  -AM  (  k  ) /XK 
Dll  55  K  I  M 
OAR  <  K ) -  0 . 0 
DO  54  I  -  i  NDA I  A 

54  OAR  (K  )  •  OAR  IKDIYII.KI  -All  <  K  >  >  **. 

55  l  i  IN  I  |  N(  II 

Dll  Ml  K  I  .  h 
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APPENDIX  C. 
(Continued) 


(MW  60  OAR<K>=0AH<K>/<XK-5  .  > 

0120  D3=o.o 

01?)  DO  68  Kd  ,M 

013?  6S  1>3=I>3+VAR  <  K  > 

0123  D3*2,#P3 

0t?4  DO  70  K-'l.M 

«1?S  70  ftlHK ,NJ)=AM(X) 

0136  80  CON  IINIJE 

0137  CALL  GDI  Ul  <  Yd  .  5  )  ,  Yd  .2)  ,-l  00  ,~1  U.?> 

0128  DO  Hi  I-l.HC-i 

01, VV  DO  8)  j=j  +  j,nc; 

0130  81  Dd.J)--fl.O 

0131  DU  82  1=5  .80-5 

013?  DO  82  J-I+l. NC 

0133  DO  8?  K--1  .  M 

0134  82  IHUDDU  .  1 )  + < AU<K . J >-AV<K . J > >#*2 

0 1 3S  DO  83  I  -'1  . NC-1 

0136  DO  83  J=)+1NC 

0137  83  Dll,  .1  >“SWRT  ( IK  1  .  T  > ) 

0138  WRIT! (6,87) 

0139  87  I-  ORMATdX  > "THE  INTERS!  I  D) STANCES " ) 

0140  DO  86  [ =1 .NC-1 

0141  DO  86  J  = I + 1 , NC 

014?  WR 1 1 E  ( 6 , 84  )  1  .  T  .  D  U  .  .7  1 

0143  84  KIRMA  I  d  X  .  "  1  ~  "  , I  1  .  2X  .  "  J=  II  ,  ?X  .  “DO  ,  J  )  =  "  ,{  1  L  ,  8 ) 

0144  86  CONTINUE 

0 1  4  S  CALL  Ct.OStdPCfO 

0146  999  ST 02 

0147  1  Nil 

014R  END* 
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